Dynamics of tachyon field with an inverse square potential in loop quantum cosmology 
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The dynamical behavior ol tachyon field with an inverse potential is investigated in loop quantum 
cosmology. It reveals that the late time behavior of tachyon field with this potential leads to a 
power-law expansion. In addition, an additional barotropic perfect fluid with the adiabatic index 
< 7 < 2 is added, and the dynamical system is shown to be an autonomous one. The stability of 
this autonomous system is discussed using phase plane analysis. There exist up to five fixed points 
with only two of them possibly stable. The two stable node (attractor) solutions are specified, 
and their cosmological indications are discussed. For the tachyon dominated solution, the further 
discussion is stretched to the possibility of considering tachyon field as a combination of two parts 
which respectively behave like dark matter and dark energy. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

The tachyon field with various potential has been stud- 
ied a lot in cosmology. Many models have been built by 
treating tachyon field as inflaton field [l|, candidate of 
dark energy [1, or a dual role of the two But 
Ref.jl] also argued that, for the class of potentials which 
V{4)) — ^ as (/) — ^ oo, radiation domination will never 
commence since the tachyon field energy density can 
at best scale as a~^. The radiation energy density would 
always redshift faster than the tachyon field. The tachyon 
field with an inverse square potential is shown to be 
able to produce a power-law expansion [5]. Coupled 
with a barotropic perfect fluid, the dynamical behavior of 
this potential has been studied in classical cosmology [_6j. 
However, the tracking solution, in which the energy den- 
sity of the field and the barotropic fiuid scales as a same 
power of a, may not be viable because its constraints on 
the adiabatic index 7. 

Our work of tachyon field cosmology is performed un- 
der the framework of loop quantum cosmology. LQC is 
a canonical quantization of homogeneous spacetime us- 
ing the techniques developed in loop quantum gravity 
(LQG). The loop quantum effects can be very well de- 
scribed by the effective theory of LQC. A modified Fried- 
mann equation is proposed and two corrections are often 
considered: the inverse volume correction and the holon- 
omy correction. However, the holonomy correction dom- 
inates over the inverse volume correction for a universe 
with a large scale factor, and thus the latter can be ne- 
glected without harm. Therefore, we only consider the 
holonomy correction in this paper. 
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Currently, tachyon matter has not been thoroughly in- 
vestigated in loop quantum cosmology (LQC). A. A. Sen 
0] generalized the description of tachyon matter in stan- 
dard cosmology to LQC under inverse volume correction. 
Xiong and Zhu 1] investigated the infiation scenario of 
a pure tachyon field with an exponential potential under 
the holonomy correction. Xiao and Zhu [8] performed 
a phenomenological analysis of tachyon warm inflation, 
in which the tachyon field with an exponential potential 
is coupled with radiation and interaction between the 
two matters were considered. However, the dynamics of 
tachyon field in LQC has not been investigated, yet. In 
this paper, we discuss the dynamics of tachyon field cou- 
pled with a barotropic perfect fluid in LQC. We focus 
on the inverse square potential, and try to explore the 
possibility of a viable dark energy model. 

The organization of this paper is as follows. In Secllll 
we try to derive the expansion law for a pure tachyonic 
matter in LQC. In Sec lIIIl we couple the fleld with a 
barotropic perfect fluid and analyze the dynamics of the 
autonomous system. The cosmological implications of 
the phase plane analysis are presented in Secs lIV Al and 
IIVBI The conclusions are made in Sec. |Vl 



II. TACHYON MATTER IN LOOP QUANTUM 
COSMOLOGY 

Based on the holonomy correction in loop quan- 
tum cosmology, the modified Friedmann equation of a 
flat(fc = 0) Friedmann-Robertson- Walker (FRW) cosmo- 
logical model is 

where H is the Hubble parameter, p and pc denote the 
matter density and critical density, respectively. We also 
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set SttG = 1 for convenience. The energy conservation 
equation is the same as the classical one, 
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where p is the pressure. Differentiate Friedmann eqaution 
with respect to time, we have 
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where 'dot' denotes the derivative with respect to the 
cosmological time t. Therefore the conditions for super- 
inflation {H > 0) are 
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where w = p/p is the equation of state. It's easy to 
see the existence of superinflation in LQC is purely an 
effect of quantum geometry, because it originates from 
the time derivative of the modification term (1 — p/ Pc)- 
The Raychaudhuri equation then becomes 
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which indicates the conditions for d > 0: 
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The regions for H > and a > are portrayed explicitly 
in FiglT] Obviously, as the matter density decreases, the 
correction term becomes less and less important, and the 
Friedmann equation as well as the conditions for a > 
and H > are consistent with classical cosmology in an 
asymptotical way. 

Now we consider the model with only tachyon field. 
Note that the tachyon field referred in this paper is just 
a scalar field with an nonquadratic kinetic term. We 
don't claim any identification with the tachyon in string 
theory. According to Sen [13] , the energy density and 
pressure of the tachyon field in a flat FRW cosmology 
can be expressed as 
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where (/> is the tachyon field, V^(0) denotes its potential, 
and the equation of state is 

V4> 



W0 



We see ranges smoothly from —1 to 0. Using the 
energy conservation Eq. ([2]), the evolution equation of 
tachyon field can be written explicitly as 




FIG. 1: a > in the region between the two solid curves. The 
region I and II correspond to _ff > 0. 



where V = dF/d^. 

Here, we are interested in the inverse square potential 



(8) 



with /? > 0. According to [5j] this potential is able to pro- 
duce a power law expansion in classical cosmology. In 
classical cosmology, one can use the Hubble parameter 
instead of the field (/) as a fundamental quantity by em- 
ploying the Hamilton- Jacobi formulation. By using Eqs. 
([T]) and ([3]) and dropping the pj Pc terms, one can obtain 
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which indicates H < 0, and therefore superinflation will 
not occur in classical cosmology. Divide both sides by (p, 
we have 



2H' 



and thus 



(10) 
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For the inverse square potential, an exact solution H ^ 
(f)^^ is found, and after some algebra one can arrive at 




where 
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The solution is inflationary when n > 1, or /? > ■y/4/3. 
Although one can have arbitrarily fast expansion with 
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an arbitrarily big n or (3, the condition for an acceler- 
ated expansion already calls for an energy scale close to a 
Plank mass Refs-jl-fllj also discussed the problem 

of tachyon model in reheating. Therefore this tachyon 
model is more suitable as a dark energy model rather 
than an inflaton model. 

However, it is difficult to find an exact solution like this 
in LQC due to the existence of the quantum modification 
term. Combining equation (1) and (3), we have 



2H 1 - p^/pc 
'3H^l-2p4,/pc 
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equation of state is P-y = — 1)P7, where 7 is a con- 
stant. Therefore, the equation of state for the whole is 

P UJ^p^ + UJjPj 

OJ — — — 

P P 

= w^n^ + oj^n-f = - 7)f]^ + (20) 

where 

Pi 1 
= — = 7 - 1, 

P7 

and Q,f, and fi^ are the fractional densities defined as 



Obviously, H > when (1 — 2p^/ p^) < 0, which means 
supcrinflation will naturally happen and it's purely an 
effect of quantum geometry as we said before. Divide 
both sides by ^, we obtain 



and therefore 
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Fortunately, we are still able to analyze its asymptotical 
behavior. When p^ ^ pc, we can neglect the derivatives 
of the correction terms since these terms, since (1— p^/pc) 
and (1 — 2p^/ Pc) will not change significantly. Then, an 
approximated power-law expansion 
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can be obtained, with 
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So we have 
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Obviously, m —> n as p^/ pc ^ 0- The quantum geome- 
try results in a different evolution of tachyon field and the 
scale factor, but the evolution will converge to the classi- 
cal one at late time when the quantum effects vanishes. 



III. WITH BAROTROPIC FLUID 

In order to obtain a viable dark energy model, we add 
a barotropic perfect fluid in our model, for which the 
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If we simply assume there is no interaction between the 
tachyon field and barotropic fluid, then their evolution 
equations are 
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and 
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Eq.Q now becomes 
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Combining Eqs. and the Friedmann equation, 

we can construct a 4-dementional autonomous system. 
To see this, we usually introduce 4 convenient variables 
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Moreover, we use iV = Ina'^ instead of the cosmological 
time t as an independent variable, therefore for any time- 
dependent function /, we have 
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Here, one should be careful that this treatment may not 
be practicable when a = 0, however, at that point, one 
can always switch back to t without causing any problem. 
With the help of new variables, the Eqs. ([T]) and ([22)1 - 
can now be expressed respectively as follows: 
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Using the above equations, differentiate x, y, z, and 
A with respect to N, the nonhnear dynamical system is 
shown in an obvious autonomous form: 
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can be fixed if the potential is specified. Note that, the 
autonomous system in LQC has one more dimension than 
the one in classical cosmology, because of the modifi- 
cation term p/pc J^- For the inverse square potential 
V{(j)) — investigated in this paper, T — 3/2. More- 

over, we find A = —ay for this potential, where 
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which means the system can be reduced into a 3- 
dimensional one: 
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Before we analyze the autonomous system described 
by Ea. ((3T|) . we should state several physical restrictions. 
First, the density should be limited and real- valued, 
therefore — 1 < a; < 1. The potential is required to be 
positive, so we have y > 0. In the analysis, we add 
the three end points of x and y for convenience, i.e. we 
set X & [—1,1], y & [0, +oo]. Moreover, the total density 
should be nonnegative and lower than the critical density, 
i.e. < z < 1. Moreover, we apply a normal restriction 
for 7, that is < 7 < 2. Therefore, w will range from —1 
to 1 due to the choice of 7. Since neither of the two parts 
should have a negative density, we have $7^,57^ G [0; !]• 
This restriction leads to (1 — z)^?/^ + x'^ < 1, and the 
equal sign is for fi^ = 1. This indicates the physically 
possible phase plane projection on the x — y plane shrinks 
as z decreases from 1 to 0. In fact, as z — >■ the phase 
plane projection becomes y^ + x^ < 1 asymptotically. 

The equilibrium points or fixed points (xe,ye, Ze, Xe) 
are solutions acquired by setting 
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According to Lyapunov's theory of stability, the stability 
of a fixed point can be determined by the property of 
the linearized system about it. The linearization is done 



by expanding Ea. ((3T|) about the fixed points and keeping 
only the linear parts. After that, one can obtain a matrix 
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at each fixed point. If the all the eigenvalues of (p3| for 
la fixed point have negative real parts, the fixed point is 
(locally) exponentially stable. However, if there exist one 
or more eigenvalues which have positive real parts, then 
the fixed point is unstable. In both stable and unstable 
cases, the fixed point is a node if the eigenvalues are all 
real, otherwise it will be a spiral point. Furthermore, the 
fixed point is called saddle point if the eigenvalues have 
both positive and negative real parts. 

Wc found up to 5 fixed points for Eg. ([31]) , their prop- 
erties are listed in Table I, where 



,0 < aj/i < 1. 



(34) 



The physical requirements mentioned before also set 
restrictions on the existence and stability of these fixed 
points. By considering 7/0;^/ y/1 — 7 < 1, we found 
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where f{j) = ^o? - 20-fa'^ + llj^o? -t- 167^^! " 7i which 
indicates the two eigenvalues of P3 with seemingly com- 
plicated square root part are definitely real and nonposi- 
tive. Note that, when 7 > c?y\, P3 doesn't exist. There 
are only four fixed points in the system, and P4, as a 
stable node, is the only attractor in the system. When 
7 = a^yl, a bifurcation occurs as the fifth fixed point, P3, 
emerges and coincides with P4. As 7 decreases from a^yf 
to 0, P4 turns into a unstable saddle and P3 becomes a 
new attractor. The position of P3 is located on a straight 
line on the x-y plane described by {x,y) = i^/J,^/J/a) 
and in the limit case 7 = 0, it will coincide with Pi. 



IV. LATE TIME EVOLUTION 

From TABLE I we can see that all the fixed points 
locate at z = 0, where the energy density vanishes and 
the LQC modification term p/pc becomes unimportant. 
We have proved in last section that the the physically 
possible phase plane projection on x — y plane becomes 
y'^ + x^ < 1 asymptotically. Therefore for the late time 
evolution, it is sufficient for us to investigate the phase 
projection on the x — y plane at small z by noting that 
z monotonously decreases as a or iV = Ina"^ increases, 
which we can easily read from Eq. ([29]). In the previous 
section, we have proved that there is at most one stable 
point in our model, then nearly all solutions will end there 
(see Figsl2]and[3]). This attracting behavior allows us to 



5 



TABLE I: Properties of fixed points 
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Coordinates 


Existence 


Eigenvalue 
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Stable node, if 7 > a''yi^; 
Unstable saddle, if 7 < a^yi^ 



investigate the properties of late time evolution near the 
attractor, regardless of the history before it. The fine- 
tuning problem can thus be waived since the same ending 



occurs for a wide range of initial conditions. Here we 
analyze two different solutions, the tachyon dominated 
solutions and tracker solutions. 
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FIG. 2: Phase plane for 7 = 0.3, a = \. The outer contour 
corresponds to y^ -\-x'^ = 1. Solutions start at z=0.01. Nearly 
all solutions end up at P3. 



A. Tracker solutions 

When -f < a^yf, P4 becomes an unstable saddle while 
^3(\/7i %/7/q^: 0) comes out as the only attractor in our 
system. The fractional density is a constant at P3, which 
means both matter scale as the same power of a. There- 
fore, we have 



lie? 



W0 



(36) 



FIG. 3: Phase plane for 7 = 4/3, a = 1. The outer contour 
corresponds to + = 1. Solutions start at z=0.01. Almost 
all solutions end up at P4. 



near the attractor, where 0o is an integration constant. 
For a same /3 or a, this solution produces a faster expan- 
sion than the one with tachyon as the only matter. It 
will lead to an eternal acceleration if 7 < min{2/3, a^yf}, 
while if 2/3 < 7 < a'^yf, a deceleration phase will def- 
initely occur at late time. Due to the choice of 7 and 
a, fl^ can range from to 1. However, Refs. 0, Q ar- 
gues that this solution cannot be a viable one, since its 
existence requires 7 < a^yf < 1. 



B. Tachyon dominated solutions 



0(t) « V7^ + 00, 



Pj (X (X a '^'^ , a{t) cx . 



(37) 



(38) 



In the case 7 > a'^yf, P4 does not exist. As the so- 
lutions converge to the only attractor P3{ayi,yi,0) in 
late time evolution, tachyon will become dominant while 
both matters are decreasing. In fact, near the attractor 
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we will have 



Pc oj ^ 0J4, a - 1, 



(j){t) « ayit, 



(39) 



P4, oc a-3(i+-*) = a-^^"y", ait) oc t^/^(^+^*) = t2/3a=,?^ 

This result is consistent with the result we derived in the 
second section. It is easy to see that if a^y^ < 2/3, the 
solutions depict an eternal acceleration. On the other 
hand, the universe will end up in deceleration if a^y^ > 
2/3. 

An interesting possibilit y is to see tachyon field as a 
combination of two parts [Tsl . which behave like a 
pressureless dust (dark matter, denoted by lower case 
DM) and a cosmological constant (dark energy, denoted 
by lower case A), respectively: 



= PDA/ + PA, 



PDM + PA, 



where 



Pdai — — ; , Pdm — u, 
1 -(^2 



and 
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In this way, dark matter and dark energy originate from 
a same scalar field, and the dynamics of tachyon field 
becomes a description of their dynamics and interaction. 
The ratio between the two parts is 



Pdm 
PA 



1 



(43) 



The proportion of dark matter rises as increases. If the 
barotropic fluid is radiation(7 = 4/3), then, by the virtue 
of the attractor solution, it is possible to have a trajectory 
that goes from the radiation dominated era to the matter 
dominated era, and then to the dark energy dominated 
regime. Radiation should dominate first, that is fi^ « 
1. Then, to have the universe dominated by matters 
described by tachyon field after the era dominated by 
radiation, that is fl^ w 1 , we need the trajectory to stay 
close to the boundary of the phase plane, i.e. (1 — z)'^y^ + 



1, or 



1 because z is very small. To have a 



sufficient long matter dominated era before dark energy 
take over, we just need the trajectory to stay close to the 
saddle point P2± after radiation's domination. Therefore 
we will rule out the trajectories which are far away from 
the boundary of the phase plane in following discussion. 

The negative/positive branch, which starts near 
P2-/P2+, will show a quite different evolution, because 
the attractor only lies in the right half of the phase plane 



and therefore the negative branch will have to get across 
to arrive at the attractor. The cosmological consequence 
will also be quite different. For the positive branch, the 
ratio of dark energy will increase as <j) decreases. If the 
attractor is inflationary (a^yf < 2/3), the allowed trajec- 
tories can go naturally into the dark energy dominated 
regime that leads to acceleration. While if the attrac- 
tor is not inflationary (a^y^ > 2/3), acceleration will not 
take place because (fi decreases almost monotonously for 
the feasible trajectories which are close to the boundary. 
On the other hand, for the negative branch, the universe 
will definitely enter a dark energy dominated regime that 
leads to acceleration when (p"^ < 2/3, or —2/3 < (j) < 2/3, 
since the trajectories will get across from the left half 
to the right half. Note that, for the negative branch, 
the universe is not expanding in the power-law way we 
described before when it first gets into the acceleration 
regime —2/3 < < 0. When = 0, dark energy domi- 
nates completely. After that, the ratio of matter increases 
again, and the final state of universe will(or will not) be 
inflationary if a^yl < 2/3(or a'^yf > 2/3). So there is 
possibility that we are just currently living in a transitory 
accelerating period and the acceleration rate can change 
according to the dynamics of the tachyon field. 



V. CONCLUSION 

Previous works on tachyon cosmology in spatially fiat 
FRW universe showed a purely tachyonic matter with 
an inverse square potential V — I3(j)^'^ leads to a power- 
law expansion In LQC scenario, although it is hard 
to find a exact solution, the expansion of universe is 
nearly a power-law one when p/pc is small. When the 
tachyon field is coupled with a barotropic perfect fluid 
with < 7 < 2, we are able to find two kinds of stable 
nodes which exist exclusively to each other and repre- 
sent different cosmological situations. The tracker solu- 
tion exists for 7 < a^yf, while the tachyon dominated 
solution exists for 7 > a'^yf. Refs.[l, Q argued that 
the tracker solution cannot be a viable one, since its ex- 
istence requires 7 < a^yf < 1. Therefore we focused 
on the tachyon dominated solution. We considered the 
tachyon field as a combination of two parts which re- 
spectively behave like dark matter and dark energy. Be- 
cause of the existence of stable node(attractor), we found 
when 7 = 4/3, it's possible to have a trajectory in which 
the universe evolves naturally from radiation dominated 
regime to matter dominated regime, and then into the 
current dark energy dominated regime through quali- 
tative discussion. The negative and positive branches, 
which respectively start from close to < and 4> > 0, 
can be interpreted into different cosmological evolutions 
and it is possible that we are just living in a transitory ac- 
celerating period, while the final stage of universe will be 
identical for the same a or /3. However, this perspective 
to see tachyon field as a combination of two parts revives 
the need of fine tuning. We also avoided the discussion of 
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the period before radiation dominated era. Further work 
can be done to bridge the gap. 
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